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Abstract 



The subject of analysis is a non-linear three-compartment model, widely used in phar- 
macological absorption studies. It has been transformed into a general form, thus leading 
automatically to an appropriate approximation. This made the absorption profile accessi- 
ble and expressions for absorption times, apparent permeabilities and equilibrium values 
were given. These findings allowed a profound analysis of results from non-linear curve 
fits and delivered the dependencies on the systems' parameters over a wide range of values. 
The results were applied to an absorption experiment with multidrug transporter-affected 
antibiotic CNV97100 on Caco-2 cell monolayers. 
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1 Introduction 



Orally administered drugs are mainly absorbed by the small intestine; they are me- 



diated upon by a variety of processes (Hunter and Hirst 1997). The drug passes 
from the lumen through the epithelial cells and the lamina propria into the blood 
stream (fig. [T] (left)). On its way it can be metabolised, transported away from the 
tract where absorption is possible or accumulate in organs other than those of treat- 
ment. 



Much experimental activity aimed at analysing the kinetic aspects of the process of 
drug absorption has been pursued recently. For better control, a variety of in-vitro 
methods of drug absorption have been developed (Balimane et al. 2000[ ). Epithelial 



^ Corresponding author. Email address: nikoOifisc .uib-csic . es, Tel: -1-34-971- 
259520, fax: +34-971-173426 



Preprint submitted to Elsevier 



12 January 2009 




V. 



apical 



cellular 



V I 

basolateral 



Fig. 1. schematic view of: intestinal wall with microvilli and a capillary embedded in the 
lamina propria (left); epithelial cells on filter (middle); simplified mathematical model with 
efflux pump indicated on apical side (right) 

cell cultures can be seeded in a mono-layer, forming the contact surface between 
two chambers (fig. [T] (middle)) and concentrations of an applied drug can be mea- 
sured over time in both chambers. Two of the well-known cell culture models are 



Caco-2 cells ( 


Artursson and Borchardt 


1997 


Artursson et al. 


cells ( 


Irvine et al. 


1999 


). 



Apart from a fundamental interest in understanding the basic mechanisms by which 
a drug is assimilated by the human body, the kinetics of drug absorption is also a 
topic of much practical interest. Detailed knowledge of this process, resulting in the 
prediction of the drug absorption profile, can be of much help in drug development 
stage ( [Eddershaw et al.[[2000[|Zhoul|2003| ). To this end, several kinetic models for 
drug absorption within the body have been established (see e.g. ( |Yu and Amidon[ 
[1999)). In this paper we analyse in detail a previously developed model, belonging 
to the category of the so-called three-compartment models ( jSkinner et al.[ 1959 



Kramer et al. , 1974) in which substances move between three volumes (e.g. the 
human organs) and active pumps are modelled by terms of non-linear fluxes. We 
will provide an analytical solution to an (adequate) approximation of the model, 
which facilitates analysis of the absorption characteristics as a whole without re- 
quiring repetitive numerical integration of the differential equations. This facilitates 
a fast and easy insight into how physiologically meaningful parameters influence 
quantities available from experiment in three-compartment models. 

The paper is organised as follows: After introducing the model and presenting its 
solution in section |2} we discuss the results and apply them to an antibiotic ab- 
sorption study carried out with the compound CNV97100 on Caco-2 cell cultures 
(section [3]). Conclusions are found in section |4j Details of the mathematical proce- 
dure are summarised in the Appendix. 



2 The model and its theoretical treatment 



Compartment models describe the behaviour of solutions or emulsions in connected 
volumes by analysing the molecule flux between them and all sources and sinks. 
When applied to drug absorption some specific simplifications must be made: Here 
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it is considered that two volumes (e.g. gastrointestinal lumen and blood plasma in- 
vivo or apical and basolateral chamber in-vitro) are connected through a third, in- 
vitro: cellular, volume. The absorbed substance should have low lipophilicity such 
that the unstirred water layer can be neglected. Furthermore it is assumed that the 
compound does not ionise and that the concentrations in the different intestinal cells 
of the mono-layer are equal (which is exact only if all cells have the same param- 
eters). With these assumptions the absorption can be seen as a transport from one 



large volume to another through a third (the cellular) volume (Gonzalez- Alvarez 



et al.[ 2005). Figure [T] (right) sketches the simplifications of the model. There is 



no spatial dependency and molecules can pass through the two cell membranes. 
The overall amount of drug molecules is considered to be constant, a hypothesis 
that assumes a closed system and that metabolism does not occur. For an in-vitro 
experiment of short duration this is a reasonable assumption. 

Passive transport across the membrane is mediated, to a first approximation, by 



the concentration gradients according to Fick's law (Fick 1855 1, which specifies a 
linear relation between the flux of particles and the concentration gradient. When 
passive absorption is accompanied by energy-consuming efflux transporters, it is 
represented by a non-linear function term in kinetic transport equations. A variety 
of transporter types could be involved in the absorption of the molecules. In our 
work we consider that the non-linear transporters are present only on one (the api- 
cal) cell membrane, but our results could be extended directly to the case that those 
transporters are located on the basolateral membrane (or even in both membranes). 
Incorporating both linear and non-linear terms, the time evolution of the amount of 
diluted molecules (Qa/c/b) in the three compartments can be described as follows: 
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(Ic) 



Qo = Qa + Qb + Qc, (Id) 



where equation (Id) stands for conservation of the overall molecule number, Qq. 
The indices denote the corresponding compartment (Apical, Cellular, Basolateral), 
Va/c/b are the respective volumes. The apical, cellular and basolateral concentra- 
tions are given respectively by <3 = Qa/Va, c = Qc/Vc and b = Qb/Vb- The passive, 
linear, diffusion terms are proportional to the concentration difference, being CIac 
and ClcB the clearances indexed with their respective membrane index. In the equa- 
tions, ^ represents the non-linear contribution due to specific efflux transporters. 
As it is an energy-consuming process, this can happen both along or against the 
gradient. 
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In the appendix we will introduce an equivalent way of writing these kinetic equa- 
tions ([T]) that will demonstrate its formal similarity with some problems in the field 
of mechanics. This mechanical analogy will invite some approximations, consist- 
ing in a linearisation of allowing us to find explicit solutions for the evolution 
of the number of molecules on each compartment. If approximated this way, any 
mass-conserving, three-compartment- model will result in a sum of three exponen- 
tials, independent from the exact form of non-linearity: 

(2a) 
(2b) 
(2c) 

where ?2 and define three time-scales and Q""^, Q"^ and Q"^ are the equilib- 
rium asymptotic quantities of the diluted substance in each compartment. These 
and the constants (A/5/C)(i/2/3) adopt different expressions, depending on the non- 
linearity. If given a specific expression for the non-linear transport terms and numer- 
ical values for the parameters, one can calculate the above constants and compare 
the result with the numerical integration of the non-linear system ([T]) or with experi- 
mental data. This explicit type of solution for a Michaelis-Menten flux J constitutes 
one of the main results of this paper and is the basis for the subsequent analysis. 
In section [3] we will carry out this program explicitly for the model and data taken 
from (Gonzalez- Alvarez et al. 2005 1, resulting in parameters as given in table|4j 
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Experimentally, it is rare to measure the complete variation of 2a (0> 2b(0 and 



Qc{t) with respect to time. A typical experiment (Lentz et al. 2000; Ruiz-Garcia 



efaLl [20021 [Faassen et aLj |2003t [Balimane et al.] |2004t [Gonzalez-Alvarez et al.j 



2005 [ ) starts by placing an initial concentration Co of a drug in either the apical or in 

the basolateral compartment. The so-called apparent permeability P"'^'' = ^^^^^ ^ 

S Co 

with Q{t) the amount of material on the receiving side, is measured in both direc- 
tions and the values are compared. 

The explicit solution Eqs.Q identifies three different characteristic time scales, ti , 
t2 and ?3, within the evolution of these concentrations. Each one of them separates 
well defined regimes in the evolution of the concentrations: if time is much smaller 
than the characteristic time scale, the corresponding exponential term comes close 
to being linear; it changes exponentially at times close to it and is almost constant at 
much larger times. This information helps the experimenter to decide if the chosen 
sampling interval is adequate or not. Furthermore, it is mathematically possible to 
observe oscillatory behaviour in the presented system, if one of the ^1,2,3 is complex. 
The condition for this to happen is shown in the appendix. 

We would like to stress that our way of approximating the problem has allowed us 
to identify these natural time scales and to find their relationship to other constants 
that are experimentally accessible. It appears that in many cases (one example in 
the next section and table [1]) one time-scale (^2) is much smaller than the inter- 
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val between measurements. Measurements of the apparent permeability are usu- 
ally carried out within a time frame of between 15-30 min and a couple of hours 
dGonzalez- Alvarez et al.] [2005| |Lentz et al.j [2000i [Yamashita eTaLj [20001 ), whereas 
t2 seems to be of the order of a few minutes. Hence, measurement times satisfy 
f » ?2 and the exponential term er'^'^ can be neglected. Analysis of the experiments 



( [Gonzalez- Alvarez et al.[ |2005[ ) indicates that transport is mediated by transporters 
with an intracellular binding site. In this case, and according to table [T] both ty and 
?3 are much larger than the measurement times, allowing us to perform the linear 
approximations e^'^'' « 1 - t/ti and e"'^'^ « I - t/^io obtain explicit expressions 
for the apparent permeability: 



pcpp 
^ BA 



1 Ml ^ A3 

SCo\ t\ ?3 



(3) 



in the case that the drug is initially delivered in the basolateral side, and 
1 



papp _ 



Bi ^ Bt, 

SCo\ ti ?3 



(4) 



when the drug is delivered in the apical side. It seems that even in the case of very 



fast- absorbing drugs (tested for permeabilities such as those in (Korjamo et al 



2007 )) times are well separated into those of around an hour and those of less than 
a minute (data not shown). However, when it occurs that the time scales are all of a 
similar order, one can easily extend eqs.[3]and|4]by the required term. 



Absorption of many drugs ( Raviv et al.[ |1990[ ) is seriously limited by P-glyco- 



protein (P-gp), the multidrug transporter. This particular protein is expressed on the 
apical membrane of intestinal epithelium cells (Tro utman and Thakker[ 2003 a|b 



Ruiz-Garcia et al.[ [2002| ). The molecule to be transported must bind with the pro- 



tein and will then be "flipped" ([ Hunter and Hirst[ 1997[ ) onto the other side of the 



membrane, where it is no longer available for the "reaction". This makes its dy- 
namics similar to enzyme reactions and is often represented by the sigmoid shape 
of a MicHAELis-MENTEN-reaction rate ( [Michaelis and Menten[[19r3l ): 



Vm determines the maximal reaction velocity, S is the surface area and Km is the 
concentration for which the velocity extends halfway towards the maximum. In 
this case of an intracellular binding site, the relevant variable is the concentration 
Qc/Vc around the binding site of the transporter, inside the cell. 

In other cases, the efflux pump possesses an extracellular binding site. Conse- 
quently, transport is determined by the drug concentration, Qa/Va, in the apical 
compartment and the corresponding Michaelis-Menten expression is: 
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J{Qa) = 



Fig. 2. Structure of antibiotic CNV97100 
SVmQaIVa 



Km + QaIVa 



(6) 



Both situations will be considered in this paper. For reasons of simplicity, we have 
considered that the efflux pumps depend on concentration on one of the two sides of 



the membrane (Ruiz-Garcia et al.l 2002 , Gonzalez- Alvarez et al. 2005 1. However. 



new results suggest that the transporter binding site for the molecule is inside the 



inner leaflet of the membrane (Hennessy and Spiers , 2007). If we consider the space 
inside the phospholipid bilayer to be an additional volume with two permeable 
walls on either side, the concentration in that volume would be between those in 
the adjacent volumes. 



In the following section, these approximate expressions will be compared with the 



experimental data from (Gonzalez- Alvarez et al. 2005). 



3 Results 



Our mathematical treatment of the evolution equations has provided us with ex- 
plicit expressions for the evolution over time of the amounts 2a/b/c(0> Eqs.(|2]),and 
the apparent permeability, Eqs.Q and (|4]). In this section we intend to determine 
the validity of our approach. Firstly, by comparing theoretical predictions with the 
full results of numerical simulations of the evolution equations, we will show that 
our approximate treatment has a range of validity that covers typical experimental 
situations. 



Next, once this validity has been established, we will use our approximation to 
make comparisons to experimental data from an antibiotic (CNV97100, see fig. [2]) 
absorption study (Gonzalez- Alvarez et al. , 2005). In this experiment a pH of 7.0 
was used in both the apical and the basolateral chamber to avoid any bias due to 
ionisation effects. If the pH were different in the two chambers it would be nec- 
essary to include these effects into the model by estimating the fraction of ionised 
and non ionised compound and to model separately the permeation of both species 
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in both directions in each chamber. The unstirred water layer had not been con- 
sidered to be the limiting diffusion step, an assumption justified by taking into 
account the molecular weight and the lipophilicity of CNV97100. This had been 



checked experimentally in situ in rats (Bermejo et al.[ 1999) and in vitro in Caco-2 



cells ( [Bermejo et al. , 2004[ ). As a consequence the three compartment model with a 
Michaelis-Menten-type flow was chosen to be an adequate picture of the underlying 
processes. 

Finally, we will determine the variation of the efflux ratio, the equilibrium concen- 
trations and the characteristic time scales for absorption with the system parame- 
ters. 



3. 1 Concentration evolution and binding site location 



We obtained expressions for the equilibrium amounts Qj^i^k^ and the constants 
(A/ B/C)(i/2/3) for non- linear transporters of the Michaelis-Menten type, in cases 
of transport mediated by both intracellular and extracellular binding. The corre- 
sponding expressions are summarised in table |4] of the Appendix as a function of 
parameters of the model. Using the numerical values of those parameters as de- 



rived in the absorption study of reference ( [Gonzalez- Alvarez et al.[ [IPOS '), we can 
extract precise numerical values for the equilibrium amounts Q^/s/c ^^'^ 
stants (A/5/C)(i/2/3). Those numerical values are also listed in table|4j Finally, the 
numerical values of the time constants ^1/2/3 are listed in table [T| Using these nu- 
merical values, extracted from a real experiment and thus corresponding to a case 
of interest, we now proceed to check the accuracy of our approximation. To this 
end, we plotted in figures [3] and 4|(for two diff'erent initial conditions) the results of 
the direct numerical integratior ^ of Eqs.([T]) and our approximation, Eqs.(|2]). The 
most noticeable feature is that, at the scale of the figures, the two approaches are 
nearly indistinguishable and, in fact, the difference is of the order of the thickness 
of the lines. We conclude that our treatment provides a simple, yet very precise, 
expression for the evolution of the amounts 2a/b/c(0 and can be used instead of 
the less transparent numerical integration of the equations. 



AB ' 



Table [T] lists the apparent permeabilities and the resulting efflux ratios, P^^a /Pa 
as obtained from Eqs.([3]|4]) in the cases of internal and external binding site and 
different initial concentrations. The analysis of the three time scales ^1/2/3, listed 
in the same table, shows that ?2 is well below the time of the first measurement 
(30min) for all cases studied and the duration of experiment ~ 2 h) satisfies 
t «; ?i,?3, hence validating the approximations that led to Eqs.(|3]j4]). The same 
analysis justifies the validity of the linear fit used in the experimental studies to 
extract the apparent permeabilities from the data. To avoid overloading the paper 



^ For this numerical integration we used a fourth-order Runge-Kutta algorithm with a time 
step of 1 s. 



7 



basolateral (donor volume) 




apical (receiver volume) 



Fig. 3. Time evolution of concentrations on either side of the cells and inside. Dotted 
Une: Numerical integration of Eq. ([T]l. Continuous line: explicit solution Eq.Q. Param- 
eters taken from ( [Gonzalez- Alvarez et al.[ 2005 1, intracellular binding site with Michaelis- 
Menten dynamics ([5]) is considered. Initial concentration Co = TSOOfxM is applied on the 
basolateral side. Right graph: first hour amplified. 

with too many results, we have omitted the time scales for basolateral loading since 
they are of similar order. 

To end the comparison with experimental data, we plot the results of the CNV97 100 
study in figure|5} superimposing on the data a line with a slope equal to the apparent 
permeability from table [T] (multiplied by 5 Co). If the antibiotic is loaded apically 
(top row in fig.|5]) the prediction of the model is good for all initial concentrations. 
If the loading is basolateral (bottom row); the prediction for low initial concen- 
tration underestimates the measured slope. As a consequence, the predicted efflux 
ratio differs from experiment at lower concentrations, which makes a more detailed 
analysis of this difference necessary. The Michaelis-Menten kinetics used seem to 
provide an insufficient representation of P-gp efflux at lower initial concentrations 
when loaded basolaterally. The model underestimates the pump's efficiency in a 
basolateral to apical set up. We stress that our approximate solution still yields very 
accurate results for the apparent permeabilities and that this observed difference is 
a direct consequence of the model or the parameters used. To make this point clear, 
we have plotted the result of the numerical integration of ([T]) for the lowest con- 
centration of SO/uM in figure |5j The deviation from the experimental data is clearly 
observable in the case of basolateral loading of the drug. Considerations of other 



efflux pathways in the P-gp transporter protein are found for example in ( Acharya 
CT1 I20061 ). 



3.2 Parameter dependence 



Once we had determined the validity of our approach, we wished to use the ex- 
plicit expressions to determine the dependence on the system parameters of some 
quantities which are of experimental interest. Here lies one of the strengths of our 
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apical (donor volume) 



cellular 



basolateral (receiver volume) 



t (h) t (h) 

Fig. 4. Same as figure [3} but initial concentration SO/jM applied apically. 

Table 1 

Time scales (when loaded apically) and apparent permeabilities predicted for different ini- 
tial concentrations and different models (internal/external binding site). 



Co 


binding 


h 




t3 




1) PZ' (cms-1) 


pcipp 
'^AB 


7500 


int 


12.3h 


6.38min 


23.6h 


6.70x10"*^ 


6.37x10-'' 


1.05 




ext 


11. 9h 


6.67min 


0.567h 


6.80x10^^ 


4.12x10-'' 


1.65 


5000 


int 


12.5h 


6.27min 


23.6h 


6.73x10"*^ 


6.26x10-'' 


1.07 




ext 


11. 9h 


6.67min 


0.567h 


6.86x10-*^ 


3.97x10-'' 


1.73 


1000 


int 


13.5h 


5.59min 


23.6h 


6.88x10"*^ 


5.57x10"'^ 


1.24 




ext 


12.5h 


6.70min 


0.567h 


7.18x10"*^ 


3.10x10"'^ 


2.32 


50 


int 


13.8h 


4.90min 


23.6h 


7.10x10"*^ 


4.88x10-'' 


1.45 




ext 


15.7h 


6.73min 


0.567h 


6.94x10 


2.07x10"'' 


3.35 



solution: In figure[6]we plot the characteristic time scale for absorption t\, the equi- 
librium concentration ratio on both cell sides beq/aeq and the efflux ratio P"^'/ /P'^b 
as a function of the clearances CIac and CIcb, the pump parameters Vm and Km and 
the initial concentration of drug Cq. Analytic formulas give access to these results 
much easier than repetitive integration throughout parameter space plus extracting 
the data from the resulting trajectories. 



Again, for reasons of simplicity, we have limited our presentation to the case of a 
secretory pump located apically with intracellular binding site; the best model ac- 
cording to analysis by ( [Gonzalez- Alvarez et al4r2005 ). As observed in this figure, 
an increase in the initial concentration Cq implies a decrease in the characteris- 
tic time ti from a finite value to a minimum value, limiting ti to a certain range. 
Raising Co increases the equilibrium concentration ratio beq/aeq. Although this ra- 
tio varies significantly, steady concentration in the basolateral site, beq, shows a 
good linear dependence with Co (not shown in the figure). Note that the efflux 
ratio P''g^ IPZb ^^^'^ decreases with increasing initial concentration, a feature sup- 
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Fig. 5. Antibiotic's concentration on the receiving side. {Top: drug loaded in apical com- 
partment, bottom: basolateral loading.) Experimental values are from CNV97100 study. 
The solid line's slope is the prediction from the theoretical solution (shown in table [T]) for 
intracellular binding, which was considered to be the model of best fit. The dashed line in 
the graphs on the right shows the numerical integration of the full system ([T]). 

ported by the experimental data, although the theoretical values deviate from the 
experimental results at low concentrations, a fact already discussed in the previous 
section. The clearance CIcb of the membrane where the pump is not situated has no 
influence on the equilibrium concentration and efflux ratios, but an increase of CIcb 
decreases the characteristic time t\, indicating a faster transport of the drug. On the 
other hand, an increase in the clearance CIac of the cell membrane where the pump 
is located has the effect of decreasing the efflux ratio and increasing the equilibrium 
concentration ratio. For large initial concentrations, Co = 1500/uM, the character- 
istic time ti shows interesting behaviour with CUc since it first increases and then 
decreases, indicating very slow drug absorption for some intermediate values of the 
clearance. 

At large concentrations, the three quantities analysed show small dependence with 
respect to the pump parameters Vm and Km, since the corresponding curves are 
almost flat. This makes it difficult to extract from the data accurate values of the 
pump parameters at those large concentrations. This suggests that lower concentra- 
tions would allow for a better experimental determination of the pump parameters - 
a practise used by experimentalists - however, we have to take into account, as dis- 
cussed above, that the accuracy of the model might worsen with decreasing concen- 
tration. In the graph, we have included negative values for Vm, which is equivalent 
to a change in the flow direction of the pump. 

Apart from these considerations, analysis of parameter dependency is the first step 
towards examining the propagation of errors into the experimentally available quan- 
tities. For example, it is clear from figure [6] (third and fourth column), that small 
differences in CIac would be nearly unnoticed, due to the rather flat curve around 
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Vj^(mol/{cm\)) (nmol/cm^') Cl^^ (10"^ cm'/s) Cl^g (10'^ cm'/s) C(,(moVcm^) 



Fig. 6. Secretory pump, intracellular binding site - Top and middle: equilibrium concen- 
tration ratio (basolateral/apical) and characteristic time (both for apical loading), bottom: 
efflux ratio J^g^ IP^j^g ■ Dependence on Vm and Km (transporter parameters) and clearances 
CIacicb- Continuous line: Co - ISOOfiM, dot-dashed line: Co = 50//M. On the very right: 
dependence on initial concentration Co- Arrows on top mark the experimentally derived 
value. (Values for the respectively fixed parameters taken from table [2]) 

its measured value (marked by the black arrow on top of the figures). On the other 
hand, a small change in CIcb yields a big variation of time scale ti. More detailed 
investigations of errors and how they bias the results have been left aside for sub- 
sequent investigations. 



4 Discussion 

Three-compartment models are widely used in drug absorption studies. The one 
treated here consists of a linear part, representing passive absorption (for example 
through cell membranes), and a non-linear part which models other means of trans- 
port, such as ABC-transporter proteins embedded into the cell membrane. In this 
work we have shown a way to transform this model into the picture of a ball in a 
well (fig. |7]), which is a general form of equations facilitating analysis of mathe- 
matical structure. The lowest point of the "well" gives the value of the equilibrium 
concentration after saturation of the process. With an adequate approximation, 
one can derive the absorption profile as the sum of three exponentials and identify 
their characteristic time scales dividing the process into phases of linear change, 
non-linear fluxes and saturation. In the phase of a near- linear profile we can provide 
explicit expressions for the apparent permeability, a quantity usually measured in 
experiments. 
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These are general results which we used on a showcase system, where the non- 
linear transport is described by a Michaelis-Menten profile - a common model of 



transporter proteins (Sharma etal.[ 2002; Volk and Schneider} 2003, Mizuarai et al 



2004j). The approach presented in this paper make efflux ratio and time scales ac- 



cessible. Both apical and basolateral drug loading can be treated by changing the 
initial conditions; we furthermore considered the possibilities of intracellular as 
well as extracellular binding sites. We analysed the importance of each physiolog- 
ical parameter in a wide range of values. 

The presented results may contribute to a better understanding of the absorption 
process and to explaining differing observations in identical experimental setups 
from a more fundamental basis. Knowledge of parameter dependencies, a funda- 
mental analysis of errors (confidence interval) and their consequences becomes 
possible. This has been left aside for later studies. The most promising experimen- 
tal setups can furthermore be predicted by treating a newly proposed model in the 
same way. 
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5 Appendix 



We will now give an overview of the used mathematical methods used and show 
the explicit solutions in the following subsections. 



The conservation law Eq. (Id) reduces the number of independent variables from 
three to two concentrations. In other words, the system described by first-order 
differential equations ([T]) only has two degrees of freedom. This implies that it can 
be replaced by a single differential equation of second order of the form: 

x = -Y{x)x + F{x). (7) 



Variable x represents a rescaled concentration in one of the compartments (the cel- 
lular compartment if depends on Qc, as when using ([5]), and the apical compart- 
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Fig. 7. Visualisation of the particle-potential picture. The position of the "ball" represents 
the concentration which will eventually end up in the minimum. 



ment if J' depends on Qa, as when using (|6])); the speed, i, and acceleration, x, are, 
respectively, the first and second derivatives of x with respect to a rescaled time s, 
and Y{x) and F{x) are functions to be described below. 



A qualitative understanding of these dynamics can be had by acknowledging that 
the previous equation corresponds to the equation of motion (Newton's second law) 
for the position x of a particle of unit mass upon which a force F{x) and a friction, 
proportional to the particle speed x act. The force F{x) and the friction coefficient 
T{x) contain all of the parameters of the system as well as the particular form of the 
non-linear flux (eqs. ([TO])). 



If we now introduce the potential function V{x) from which the force derives as 
dV 

F(x) = , the evolution of x can be visualised as the relaxation of a ball rolling 

dx 

downwards within a well of shape Vix) under the combined effects of gravity and 
friction. It is known from mechanics that the particle will eventually stop at the 
minimum of the potential (stable equilibrium state) and that relaxation towards this 
final state will proceed via damped oscillations or monotonously, depending on 
the relative strength of the friction and potential contributions. Figure |7] visualises 
this ball-in-a-well approach. For very general non-linear transporters, including the 
Michaelis-Menten form, Eqs.(5j6), used in this paper, the corresponding potential 
function V(x) displays a single minimum, although its exact shape depends both on 
the linear and the non-linear terms in Eqs. ([T]). 



Once this mechanical simile has been realised, an approximation appears quite nat- 
ural: the potential V(x) is replaced by the parabolic approximation around its mini- 
mum, equilibrium, value x^q, i.e. approximate V(x) ~ V{Xeq)+\{x-XeqY . Similarly, 
the friction coefficient is replaced by its value at this minimum Y{x) ~ T{Xeq) = F^^. 
With these approximations, Newton's equation (|7]) becomes the linear equation 
X = -TeqX - o/{x - Xeg), wMch dcscribcs the damped pendulum of frequency co. 
The solution of this equation, as found in many elementary books of mechanics, 
can be written as the sum of exponential functions of time. One can then undo the 
change of variables and write the solution in full as done above (eqs.|2]) 
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5. 1 Parameters and simplifications for three-compartment system 



The conservation law (fTdl) allows the elimination of one of the three equations of the 



set ( la Ic I. We have chosen to eliminate either Qa or Qc, depending on whether the 



non-linearity depends on Qc or Qa, respectively. More precisely: if ^ = J'iQc), 

Qcit) 

we define the dimensionless normalised concentration x(t) = where A/" is a 

VcN 

normalisation constant with units of concentration to be specified later; in the other 

case, when = J'{Qa), we define x{t) = in both cases, we also define 

VaN 

y(t) = YT^- The initial drug concentration Co = Qo/Vi (being V[ = Va or Vj = Vb 
VbN 

according to whether the drug is initially loaded on the apical or the basolateral 
volume, respectively) is also rescaled to cq = — . We finally define a rescaled 

CI 

dimensionless time variable s = 1. It turns out that, with these definitions, the 

resulting two-degree system can be written in the common form: 

dx(s) 

x(s) = — — = aux + any + an + j(x) (8a) 
ds 

dy(s) 

y(s) = — — = a2ix + a22y + «23- (8b) 
ds 

Vb Vb 
with j = —z—-—JiCly{) in one case, and j = - .r ^^Qc) in the other. The 

CIcbVaN CIcbVcN 
dimensionless constants a,y depend on the clearances, volumes, overall concentra- 
tion, as specified in table [3j 



By diff"erentiating x{s) with respect to s using Eq.(8a) and, in the resulting expres 



sion, replacing y(s) by Eq.([8b|) and y{s) by its isolation from Eq.([8a|), we get the 
form 

X = -r (x) X + F(x), (9) 



with friction coefficient Y{x) and force F{x) given by: 

r{x) = To - f(x), F{x) = a -fix - a22j(.x). (10) 



Fq, a and fi are dimensionless, positively-defined constants, whose relation to the 
coefficients aij are also detailed in table [3] 
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5.2 Coefficients in case o/Michaelis-Menten efflux transporter 



This presentation has so far been very general. In the case of the Michaelis-Menten 
form for the non-linear efflux transport function ^, formulas (|5]) or (|6]), we identify 
Km as a characteristic concentration and simply adopt the normalisation constant 
N = Km- Therefore the current is 

jXx) = y-^ (11) 
I + X 



where y is written up in table [3} The friction coefflcient and force are: 

y X 

r(x) = ro-- -T, F{x) = a - fix - ya22- • (12) 

(I + xY i + X 



Note that the initial concentration Co is contained only in a and that the influence 
of the efflux transporter is found completely in y. Yq, a and /3 are independent of 
the location of the pump. The values used for numerical calculations are taken from 
( [Gonzalez- Alvarez et al.[ |2005| ) and can be seen in table [2| 



The potential from which the force F{x) = -^^j^ derives is: 

V(x) = (-a + ya22)x +2"^ ~ '^^22 ln(l -1- x) . (13) 



The equilibrium concentration Xeq is the minimum of the potential, obtained by set- 
ting the force in Eq.([T2]) equal to zero. In the present case of a Michaelis-Menten 
type non-linearity, the solution is obtained without using any further approxima- 
tions or simplifications other than stated for the original model: 



a - P - 7^22 + ^j^ctfi + {a - /3 - ya^iY 



Xeq = ' . (14) 



Inserted into ( 8b I the equilibrium value for y(s) is the solution of y(s) = 0: 



yeq = (15) 

«22 



Through the conservation law one gets Zeq 
Go K Vb 



•'eq 



- Xeq— - yeq— , (16) 
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Fig. 8. Continuous line: Original potential Eq. ( 13 1, dashed line: approximated parabola. 
For different amounts of loaded drug. Differences are of the order of the line width. 

where, to unify notation, we have labelled the volumes' meaning = Vc and 
= Va for intracellular binding, and = Va and = Vc for extracellular 
binding. 



5.3 The linearised model 



The potential V{x) , Eq. ( 13 1, has a single minimum and it can be approximated by 

a parabola around this minimum V{x) = Veq + ^ V"^ {x - . The second derivative 
with respect to x is: 



0-221 



(17) 



The approximated and original potentials are drawn in figure [8]in the case of apical 
loading of a system with intracellular binding site. From this figure, it can be seen 
that the parabolic form constitutes an excellent approximation for a wide range of 
initial concentrations Cq. 



Furthermore, we also approximate the friction coefficient Y(x) by its value at equi- 

The resulting hnear differential equation x = -TeqX - 



librium F, 



eq 



Fn - 



(l+JTf,) 



CO (x - Xeq) has the solution (Boyce and DiPrima 2001 ): 

x{s) = Cie-''^' + C2e-''^^ 



+ X 



eq 



(18) 



with time constants: 
2 



Tl 



and T2 



^eq + -^r^ 



(19) 



The coefficients Ci and C2 are determined by the initial conditions Xq and io, the 



latter being determined through Eq.( 8b ) by Xq and as io = 0.11X0+ ai2yo+ an 
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The resulting formulas are summarised in table |5j The initial values xq and yo de- 
pend on the particular model used. For example, for intracellular binding we have 
^0 = = and vn = %^^. This latter value, in turn, has to be adjusted ac- 

VcKm " VbKm ■> 

cording to where the drug is loaded: for apical loading yo = 0, while for basolateral 
loading yo = Note that ti and T2 will be complex if Yeq < 2oj. In this case the 
system would relax to the steady state by performing damped oscillations. How- 
ever, for the values of the parameters drawn from experiments, this case does not 
arise, ri and T2 are positive real numbers and the decay to the equilibrium state is 
governed by real exponentials. 



The time evolution of y can be obtained by a direct integration of the linear equation 



(8b): 



y(s) = e' 



yo + 



^ ds'e-"''-'' {a2ix{s') + a23) 



(20) 



which yields 



(21) 



with a new time constant: 

T3 = -— (22) 
«22 



and coefficients ^1/2/3, whose relation to other constants is detailed in tablejs} The 
evolution of z{s) is obtained by means of the conservation law. Finally, when we 
undo the changes of variables we obtain coefficients as given in table|4]and together 
with Eq.( 19 1 we derive Eqs. ([2]). It is worth recalling that the time scales in real units 
are 



Vn 



CI 



Vn 



CB 



CI 



Vn 



CB 



CI 



-T3 



(23) 



CB 
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Table 2 

Experimental parameters drawn from ( Gonzalez- Alvarez et al. 2005 1 used for the calcula- 
tions. 



Parameter 


Measured Value 


CIac 


14.49 X 10-^ cm^/s 


ClcB 


3.528 X 10-^cmV^ 


Vm 


6.17 X lO'^^mol/icm^s) 


Km 


0.376mo//cm^ 


S 


4.2cm^ 


Va 


Icnr" 


Vb 


Scrrr" 


Vc 


0.073 Scm^ 



Table 3 



Coefficients for equations ^ using parameter values from (Gonzalez- Alvarez et al. 2005 1 



Parameter 


INTRAcellular 


numerical 


EXTRAcellular 


numerical 


an 
an 
au 


\Vb , CIac ( Vb , Vs\l 
[Vc ClcB \Vc Va)\ 

Vb a CIac Vb\ 
VcV CIcbVa) 

CIac VbV, ^_ 

ClcB VaVc'-O 


-213.8 
-209.8 
2.220xl0^y/C() 


CIac (Vb ,Vb\ 
ClcB \Vc Va) 

CIac V'b 

ClcB VaVc 

CIac VbV, 
CIcB VaVc^O 


-173.1 
-250.4 
2.220xlO*^y/Co 


ail 


1 


1 


Va 
Vc 


-27.10 


ail 


-1 


-1 




-41.65 


023 








Vi 


3.604xlO*^V/Co 


7 


SV„ Vb 
CIcbKm Vc 


-79.41 


SVm Vb 
CIcbKm Va 


2.930 



To 


-an - a22 - 


. Vb , CIac ( Vb . Vb\ , , 
- Vc ClcB \Vc Va)^ ^ 


214.8 


a 


ai2<323 - fl^is^Jia - 


_ CIac VbV, 
" CIcbVaVc^'^ 


2.220xl0^y/Co 


P 


^11^22 ~ ^12^21 - 


- ?iciv?vc(y^^yB^vc) 


423.6 
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0.348 
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Ci 

Di 

D2 
D3 



(xo - Xeq + r2 (ai 



-eq + T2 [anXQ + ai23'o + «13 + 



Xo 



^eq 



C, 



Jo 



aiiTi + 1 

a2lC2T2 
a22T2 + 1 
■ yeq -D1-D2 
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